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Abstract. Aim of this work is the calibration of a SPECT/CT scanner (SYMBIA T2 - SIEMENS) for activity quantification in 
" m Tc-macroaggregated albumin ( 99m Tc-MAA) clinical acquisitions used for radiation dosimetry before radioembolization. 
Materials and methods: SPECT/CT calibration was performed using the methodology proposed by Zeintl et al [1]. A 
calibration factor converting cps (counts per second) to activity was determined using a uniform radioactive solution in a 
cylindrical water phantom. Recovery coefficients, background variability and contrast were obtained from acquisitions of an IEC 
NEMA torso-phantom with hot sphere inserts. " m Tc was used to prepare the radioactive solutions and different sphere-to- 
background ratio (Rsb). The phantom was acquired using the 99m Tc-MAA clinical protocol. Images were reconstructed with the 
proprietary OSEM 3D algorithm (FLASH 3D, Siemens Healthcare), scatter and attenuation correction. The phantom results were 
validated on 99m Tc-MAA in vivo patient acquisitions. 

Results: The calibration factor for " m Tc was 464.8 cps/MBq. Recovery coefficients vary with sphere volumes from 90% to 
20%. Hot sphere contrast and activity quantification improve with higher iterations and subsets while the sphere to background 
ratio and the background variability increase. A balance between a satisfactory activity recovery and a low background 
variability (<12%) was obtained with 8 iterations, 8 subsets and no filtering. Hot spheres activity was determined with an 
accuracy <7%. Activity in liver from 99m Tc-MAA patient acquisitions was 8% (mean value) which was less than the 
administrated activity to the patients. This underestimation was attributed to some extrahepatic shunt. 

Conclusions: Activity recovery using the tested scanner can be achieved within an accuracy of 7%. With a simple scanner 
calibration procedure, the parameters of the clinical " m Tc protocols can be optimized improving image quality and activity 
recovery according to the clinical requirements. 
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1. INTRODUCTION 

S PECT has traditionally considered a non-quantitative 
imaging modality. However, the latest generation of 
SPECT/CT scanners has made quantitative SPECT 
practicable thanks to the implementation of new 
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computational features such as attenuation, scatter, 
resolution recovery and dead time corrections. Image 
reconstruction, which was traditionally performed by 
filtered back projection, is now based on iterative methods. 
OSEM (Ordered Subset Expectation Maximization) 
algorithm is the method of choice used in hybrid scanner 
technology. In iterative methods image reconstruction is 
determined by different parameters such as iteration 
numbers, subsets and post-smoothing filter. As 
demonstrated in previous studies [2,3,4] the combination of 
iteration/subsets and post-smoothing affect the activity 
quantification and the image noise. Increasing the iterations 
and subsets improves the activity recovery but make images 
noisy and not applicable in clinical practice. Therefore, an 
appropriate knowledge of the quantification capability of 
the scanner used in daily clinical practice is mandatory to 
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assess a correct activity quantification and the highest 
image quality. As shown in literature [5], quantitative 
SPECT may have various clinical applications. In vivo 
quantification of " m Tc acquisitions can be achieved with an 
accuracy of about ±10% as reported in [5] and references 
therein; more challenging is quantitative SPECT with other 
radioisotopes. From our side we used SPECT quantification 
applied to 99m Tc-MAA acquisition to estimate radiation 
dosimetry in selective internal radionuclide therapy before 
radioembolization. For this reason, we investigated the 
quantification performance of our SPECT/CT scanner 
(SYMBIA T2 - SIEMENS). The aim of the work was the 
scanner calibration and the absolute activity quantification 
(with 99m Tc) both in phantom and in vivo acquistions. To 
perform the study, we used the methodology proposed by 
Zeintl et al [1]. The procedure consists of the calibration of 
the scanner which was achieved by phantom studies. 
Accuracy in activity definition is assessed using appropriate 
torso-phantom acquisitions with hot spheres simulating 
lesions of different dimensions. Phantom results were used 
to evaluate the activity in liver lesions for 8 99m Tc-MAA 
patient acquisitions. 

2. Materials and Methods 

2.7. Methodology 

We used a uniform cylindrical phantom with a volume of 
5680 mL and an IEC NEMA torso-phantom [6] containing 
six spheres of different volumes (25.5 mL, 11.9 mL, 5.5 
mL, 2.5 mL, 1.5 mL, 0.5 mL) and a lung insert. " m Tc was 
used to prepare the radioactive solutions filling the 
phantoms and the spheres. Phantom acquisitions were 
performed on a SPECT/CT SYMBIA T2 scanner 
(SIEMENS Healthcare) equipped with a low-energy and 
high-resolution collimator. Images were reconstructed using 
the proprietary OSEM 3D algorithm with resolution 
recovery (FLASH 3D, Siemens Healthcare), CT-based 
attenuation correction, energy window-based scatter. The 
attenuation map was generated from a CT scan of the 
phantoms with 130 kV, 30 mAs and a reconstruction kernel 
B08s (Siemens Healthcare). Phantom volume analysis and 
count statistics were performed using the proprietary 
software Volumetric Analysis (Siemens Healthcare). 

The work is subdivided in four steps: 

1) the SPECT/CT calibration using the uniform cylindrical 
phantom; 

2) the activity quantification using the IEC phantom with 
hot spheres and hot background; 

3) the activity quantification using the IEC phantom with a 
hot insert consisting of a 10 mL vial containing 9.99 MBq 
of " m Tc and an orange fruit injected with 1 mL of " m Tc 
(15.9 Mbq). For simplicity in this test we used a cold 
background. The radioactive orange fruit was inserted with 
the aim to test a non-uniform activity distribution; 

4) the phantom results were extended to 99m Tc-MAA patient 
acquisitions for quantification of the activity accumulated 
in liver lesions. 
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2.2. SPECT/CT calibration 

System calibration consists in the definition of an 
experimental factor ( S ) relating the count rate (cps) and the 
activity (MBq) contained in a tridimensional volume. This 
can be achieved by acquisition of a phantom filled with a 
known activity concentration. In our case, system 
calibration was performed using a cylindric phantom of 
5680 mL filled with a uniform dilution of " m Tc (196.1 
MBq, about 127 million of total counts). An isocontour 
acquisition was performed with 120 projections (10 sec per 
projection), 128x128 matrix, 2 separate energy windows for 
the acquisition of the photopeak and the lower scatter range 
(amplitude 20%). SPECT data were reconstructed with 
OSEM-3D algorithm, resolution recovery (FLASH 3D), 8 
iterations, 5 subsets, 8.4 mm Gaussian filter, scatter and 
attenuation correction. These are the standard parameters 
used for clinical data reconstruction. 

The calibration factor S of the system is calculated by the 
formula: 



Where Ct are the total counts per unit time (cps) in the 
phantom and At is the actual activity at the calibration time. 
Counts and activity are corrected for radioactive decay. 
Total counts were determined drawing a VOI around the 
phantom boundaries. 

2.3. Activity quantification 

To evaluate the accuracy of SPECT quantification we 
performed a phantom study using the IEC NEMA phantom 
with six spheres of different volumes (from 0.5 mL to 25.5 
mL). Different sphere-to-background ratio (R s /b) were used 
(61:1, 44:1, 33:1, 10:1). The data were acquired with the 
phantom centred in the machine field of view using an 
isocontour acquisition over 360 degrees with 120 
projections (10 sec each). We used a LEHR collimator with 
4.8 mm pixel. SPECT data were reconstructed using a 
OSEM-3D algorithm as specified before. Four different 
reconstructions were performed with the following 
iterations/subsets/filter settings: 8it 5subsets 8.4 mm filter, 
8it 5subsets No filter, 8it 8subsets No filter, 8it 16subsets 
No filter. The data was corrected for scatter and attenuation. 
The counts in the phantom spheres were calculated drawing 
a spherical VOI superimposed to the sphere volume 
determined following the CT boundaries. The absolute 
activity Aj for an object of size Sj was calculated from the 
relation: 

*1 = C, /s «> 

Where Cj is the counting rate in the object j and S is the 
calibration factor calculated in the previous section. The 
sphere-to-background ratio (R s /b) and the background 
variability (CV) were calculated as: 

Rs/b = Cmax /background (mean value) (2) 

CV— SD/background (mean value) (%) (3) 

Cmax is the maximum voxel value in the sphere of 25.5 mL. 
The background mean value is determined as the mean 
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value from 10 ROIs placed in the background and the SD is 
the mean SD of these 10 ROIs. 

The hot sphere contrast is defined as: 

( c sph/ c bk)—l 


C S ph ~ 1 


(4) 


(Csph/Cbk) 1 

c sp h is the measured activity concentration in the sphere 
with maximum diameter, Cbk is the background activity 
concentration calculated as the mean value of 10 ROIs 
placed in the background. C sph and C bk are the true 
concentration in the spheres and in the background. 

To quantify the spillover of activity outside the sphere 
boundaries and the consequent activity underestimation, we 
determined the recovery coefficients as: 


RC 


A sph 
Atrue 


(5) 


Where A sp h is the activity measured in the spheres and A tru e 
is the true activity filling the spheres. 


2.4. In vivo patient analysis 

The phantom results were extended to 99m Tc-MAA patient 
acquisitions. MAA examinations were performed prior to 
90-Y radioembolization with an injection of " m Tc in the 
range (150-^180) MBq for each liver branch. 99m Tc-MAA 
was injected via intra-arterial administration in one or both 
liver lobes depending on disease extension. In absence of 
extra-hepatic shunts activity injected in the liver can be 
measured from counts collected in the tumour region and 
applying Equation (1). Lesions were contoured drawing a 
VOI as shown in Figure 1 and the accumulated counts were 
determined using the Volumetric Analysis Software 
(Siemens Healthcare). The contouring algorithm is based on 
an adaptive thresholding method based on the source to 
background ratio. Contouring thresholds were calculated as 
a percentage of the maximum voxel value. Patient data was 
acquired using 120 projections (360-degree acquisition 
range), two energy windows for photo-peak and lower 
scatter acquisition, attenuation map from CT scan. SPECT 
field of view was centred on the patient abdominal region 
containing the liver. Image data was reconstructed with 
OSEM-3D algorithm (Flash 3D, Siemens Healthcare) with 
8 iterations/8subsets, no Gaussian filter, CT-bases 
attenuation correction and energy window-based scatter 
correction. 



Figure 1. 99m Tc-MAA acquisitions. 


3. Results 

3.1. SPECT/CT calibration 

The calibration factor obtained from the uniform cylindrical 
phantom was S: 464.8 cps/MBq. This value was used for 
absolute activity calculation applying Equation (1). 

3.2. Activity quantification (with hot spheres and hot 
background) 

The effects of object size and reconstruction parameters are 
shown in Figure 2 and Figure 3. Figure 2 reports the 
recovery coefficients for different sphere volumes, 
R s/b =10.5 and four combinations of iterations/subsets/filter. 
Recovery coefficients vary from 90% for volumes > 25 mL 
to 20% for volumes up to 0.5 mL. Activity recovery 
improved if no smoothing filter and a higher number of 
updates (iterations x subsets) are applied in the 
reconstruction process. 


-*-8i5s 8.4mm filter -#-8i 5s NO filter — 8i8s NO filter — 8il6s NO filter 



Figure 2. Recovery coefficient for a sphere-to-background 
ratio R s /b =10.5 and different iterations, subsets and 
Gaussian filter. 

Figure 3 and Table 1 report the recovery coefficients for 
different R s /b and 8iterations/8subsets/No filter. 
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Figure 3. Recovery coefficient for different ratio R s /b (8i 
8subsets and No Filter). 
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Sphere 

volumes 

(mL) 

Rs/b 

10.6 

32.8 

44.4 

61.1 

Mean 

values 

25.5 

84.2% 

±4.1% 

81.7% 

±4.1% 

85.7% 

±4.1% 

80.4% 

±4.1% 

83.0%±4.1% 

11.6 

77.6% 

±4.1% 

72.1% 

± 4.2% 

77.7% 

± 4.2% 

74.6% 

± 4.2% 

75.5%±4.2% 

5.5 

73.7% 

± 4.3% 

69.9% 

± 4.4% 

66.1% 

± 4.3% 

69.4% 

± 4.4% 

69.8%±4.3% 

2.55 

52.4% 

± 4.4% 

61.3% 

± 4.7% 

62.9% 

±4.1% 

56.2% 

± 4.7% 

58.2%±4.6% 

1.15 

35.3% 

± 4.5% 

50.1% 

±5.1% 

49.2% 

±4.1% 

42.2% 

±5.1% 

44.2%±5.1% 

0.5 

21.5% 

± 4.4% 

26.1% 

± 5.4% 

36.6% 

±4.1% 

24.5% 

±5.5% 

27.1%±5.3% 


Table 1. Recovery Coefficients for 8it/8subsets/ no filter 
and different R s /b. 

The Sphere-to-background ratio, the background variability 
and the hot sphere contrast are shown in Figures 4-6. 
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Figure 4. Sphere-to-background ratio for different 
reconstruction parameter sets. 
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Figure 5. Background variability for different ratio R and 
reconstruction parameter sets. 

Activity quantification and hot sphere contrast improve for 
higher update number (iteration number x subset number) 
and no Gaussian filter included in the image reconstruction 
(see Figure 3 and Figure 6). Better quantification is possible 
for larger sphere volumes (> 25 mL) where the effect of 
activity spill in and spill out are less important. The ratio 
R s /b is less influenced by the update numbers while the 
background variability increases with the iteration/subset 


number (see Figure 5). A balance between a satisfactory 
activity recovery and a low background variability (<12%) 
was obtained with 8 iterations, 8 subsets and no filtering. 
Applying these reconstruction parameters, the recovered 
activities in the hot spheres were determined for all 
phantom acquisition with a mean error <7.2%. From this 
analysis, we conclude that the optimal reconstruction 
parameters are: OSEM-3D algorithm, resolution recovery 
64 update number (8iterations and 8 subsets) and no post¬ 
smoothing filtering. 


R_SB=61.1 R_SB=44.4 R_SB=32.8 R_SB=10.5 



Figure 6. Hot sphere contrast. 

3.3. Activity quantification (with vial hot insert and 
radioactive orange fruit) 

The activity measured in the vial was 9.62 MBq with a 
difference of -3% with respect to the injected activity. The 
activity in the orange fruit was 14.8 MBq (7% less than the 
injected activity). The activity in the orange fruit was 
slightly inferior to the injected activity because of some 
leakage from the injection site. 

3.4. Patient data analysis 

The recovered activity of 8 in vivo patient 99m Tc-MAA 
acquisitions are reported in Table 2. Patient images were 
reconstructed using the optimized parameter settings 
obtained with the previous phantom studies (8iterations, 
8subsets and no filter). The liver lesions were contoured as 
shown in Figure 1 and the total counts were recorded. 
Absolute activity was then calculated from equation (1) and 
corrected with the recovery coefficients extracted from 
Table 1. 


Patients 

Total counts in 

the liver 

Measured 

activity 

(MBq) 

Injected 

activity (MBq) 

Difference (%) 

Patient 1 

148130656 

320 

350 

-9% 

Patient 2 

160030464 

346 

380 

-9% 

Patient 3 

97476544 

211 

230 

-8% 

Patient 4 

71863680 

155 

170 

-9% 

Patient 5 

57551472 

124 

140 

-11% 

Patient 6 

90944672 

197 

210 

-6% 

Patient 7 

100781472 

218 

230 

-5% 

Patient 8 

64042276 

138 

150 

-8% 


Table 2. Injected and Measured activity in 99m Tc-MAA 
liver acquisitions. 
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The injected activity and the residual activity in the syringe 
were measured respectively before and after injection in a 
dose calibrator (±1%). Residual activity in the intra-arterial 
catheter was also evaluated. The difference between 
injected and measured activity was in the range (-11% A - 
5%). This discrepancy can be explained by the fact that a 
small part of macro-aggregates tends to separate from 
99m Tc, therefore a percentage of activity flow from the liver 
towards other abdominal organs causing a loss of activity in 
the liver lesions. 

4. Discussion and Conclusions 

In this work we tested the quantitative performance of our 
SPECT/CT scanner using the procedure suggested in [1]. 
Our results show that activity quantification in " m Tc 
phantom acquisitions with hot inserts can be achieved with 
an accuracy of about 7%. This value is comparable with the 
accuracy measured by Zeintl et al. in [1]. Since activity 
quantification is affected by partial volume effects, the 
recovery coefficients for different source to background 
ratios were determined with an accuracy within 5%. The 
hot sphere phantom was acquired using the 99m Tc-MAA 
clinical protocol. Optimized reconstruction parameters were 
determined (8iterations/8subsets/no filter) and then applied 
to reconstruction of 99m Tc-MAA in vivo liver acquisitions. 
The recovered activity measured in the patient livers 
resulted in an underestimation of about -8% (averaged over 
all patients) with respect to the actual injected activity. We 
attributed this underestimation to some activity flow 
towards the abdominal regions outside the SPECT field of 
view. The results obtained with patient acquisitions were in 
accordance with those presented in a similar work by 
Gnesin S. et al in [7]. In summary with this study we 
verified that our scanner may quantify " m Tc SPECT/CT 
acquisitions with a good accuracy using a simple calibration 
method. Here we used the spect calibration to optimize the 
clinical protocol of 99m Tc-MAA in-vivo liver acquisition 
and to quantify the activity in liver lesions: for the future, 
the same method might be extended to optimize other 
clinical " m Tc protocols. 
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